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ABSTRACT 

We present a program designed to obtain age-rotation measurements of solar- 
type dwarfs to be used in the calibration of gyrochronology relations at ages of 
several Gyr. This is a region of parameter space crucial for the large-scale study 
of the Milky Way, and where the only constraint available today is that provided 
by the Sun. Our program takes advantage of a set of wide binaries selected so 
that one component is an evolved star and the other is a main-sequence star of 
FGK type. In this way, we obtain the age of the system from the evolved star, 
while the rotational properties of the main sequence component provides the 
information relevant for gyrochronology regarding the spin-down of solar-type 
stars. By mining currently available catalogs of wide binaries, we assemble a 
sample of 37 pairs well positioned for our purposes: 19 with turnoff or subgiant 
primaries, and 18 with white dwarf components. Using high-resolution optical 
spectroscopy, we measure precise stellar parameters for a subset of 15 of the pairs 
with turnoff /subgiant components, and use these to derive isochronal ages for the 
corresponding systems. Ages for 16 of the 18 pairs with white dwarf components 
are taken from the literature. The ages of this initial sample of 31 wide binaries 
range from 1 to 9 Gyr, with precisions better than ~ 20% for almost half of 
these systems. When combined with measurements of the rotation period of 
their main sequence components, these wide binary systems would potentially 
provide a similar number of points useful for the calibration of gyrochronology 
relations at very old ages. 

Subject headings: Galaxy: stellar content — stars: fundamental parameters — stars: 
binaries: general — stars: evolution — stars: rotation 
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Introduction 



Understanding the detailed formation and evolution of systems such as star clusters 
and galaxies requires knowledge of the ages of their constituent stars with sufficient 
precision. The derived star formation histories of the populations of stars in all types of 
galaxie s provide crucial constraint s on models of the formation of structure on cosm ological 



scales ( ISpringel fc Hernquist 



2003 



Nagamine et al. 



2005 



Conroy fc Wechsler 



20091 ). In the 



context of our own Milky Way, for example, the various theoretical models for the formation 
and evolution of the disk of t he Galaxy make di f ferent predictions f or the age distributio ns 



of th in- and thick-disk stars (lAbadi et al. 



2003 



Brook et al 



2007 



Schonrich fc Binney 



20091 ). Even though large amounts of information on the chemistry and kinematics of these 
objects are and will continue to become available, without a proper knowledge of the ages 
of thin- and thick-disk members it will not be possible to determine which are the more 
realistic scenarios. Similarly, age-dating of the numerous stellar streams being discovered 
by modern surveys would help trace back their origin and therefore contribute to a better 
understanding of the recent accretion history of the Galaxy. Exoplanet research could also 
benefit from precise determinations of stellar ages since the ages of host stars can be used as 
constrains for dynamical models of planetary systems and migration. Moreover, the ages of 
host sta rs would be cruci al for investigations of biological evolution in potentially habitable 



planets (ISoderblom 



2010h . 



Unfortunately, the determination of stellar ages is possible only in a few cases, typically 
for sets of stars that can be safely assumed to be more or less coeval (e.g., clusters, 
associations, etc.) and for very nearby stars for which a wealth of accurate information 
such as distances, luminosities, and detailed chemical compositions are also available. For 
plain and simple isolated stars in the Galactic field just outside the very immediate solar 
neighborhood, reliable ages are simply not possible yet. This lack of knowledge of the 
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age distribution of stars in the Milky Way thus constitutes a severe hmitation under the 
perspective of cosmology and galaxy formation, as it prevents an adequate investigation 
of the details of the formation history of our Galaxy. Upcoming experiments such as Gaia 
and LSST aim to change this situation by measuring the properties of many millions of 
stars with very high accuracy but, as of today, stellar ages will still need to be derived 
from isochrone fitting. While historically this has been the most successful method for 
deriving stellar ages in many branches of astronomy, it suffers from limitations that become 
particularly restrictive when dealing with the needs of a more ambitious mapping of 
the star formation histories of the components of the Milky Way. In this context, some 
of the limitations of iso chronal ages include their high sensitivity to errors in distance 



( ISoderblom et al. 



of isochrone points 



2005 



da Silva et al. 



20051) and stellar parameters, systematic biases arising from the samp ling 



Pont fc Ever 



2004 



Nordstrom et al. 



2004 



J0rgensen fc Lindegren 



20061), the impact of poorly understood stellar evolutionary processes 



such as microscopic diffusion, convective overshoot, gravitational se ttling, and others 



( IVandenBerg et al. 



2002 



Demarque et al. 



2004 



Michaud et al. 



20101 ). and the unavoidable 



fact that the method is most sensitive to evolved, thus relatively massive, stars, and not so 
much to stars on the main sequence (MS). Thus the importance of seriously exploring any 
alternative to the use of stellar models for deriving ages of field stars. 



One o f such alternatives is that offered by gyro chronology ( ISkumanichl Il972l : iBarnes 



2003 



20071 ). The method takes advantage of the fact that MS stars of FGK-types are known 



rates of stars in a sequence of clusters of dii 


(— ^ 

'erent ages (iStauffer et al 


. 1987 




Radick et al. 


1995 


1 


Irwin et al. 


200 


6; 


Messina et al. 


2008; 


Me] 


3om et al. 


2009 




Coll 


ier Cameron et al. 




2009 


Hartman et al. 


2009 




2010 




Delorme et al. 


2011 


1 


Meibom et al. 


2011a 


). Thus, if it 



were possible to map and quantitatively calibrate with precision such spin down from very 
young to very old ages, it should in principle be possible to apply those gyrochronology 
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relations to single stars in the field and, by measuring their surface rotation rates and colors 
(a proxy for their masses), readily obtain their ages. 

While there exists a growing body of observational evidence in support of the possibility 
of gyro chronology, a solid theoretical understanding still needs to be fully developed. The 
physics at the heart of the method involves the self-regulated interaction between stellar 
phenomena that are not easily modeled, such as differential internal rotation, magnetic fields 
interior to stars, and magnetized stellar winds that are able to carry angular momentum 
away. For a discussion of these pro cesses and the i r inte raction in the context of a theoretical 



framework for gyrochronology, see 



Barnes! (12003 



2OIOI ). In the present paper, we approach 



the subject from the perspective of the empirical calibration of gyrochronology relations for 
use as a tool for the age dating of field stars. 

Until very recently, reliably calibrated gyrochronology relations existed only for stars 
younger than about 0.5 Gyr, for the simple reason that such age corresponded to the 
oldest clusters for which rotation periods had been measured for a large enough number 
of member stars. The only constraint beyond this point corresponded to the Sun, thus 
preventing a good calibration of gyrochronology relations for stars a few Gyr old. Although 
there was never a lack of older clusters to extend these relations to correspondingly older 
ages, the limitation was of a technical nature as older stars become progressively less active, 
thus complicating the task of measuring the already small photometric modulation due to 
starspots coming into and out of the line-of-sight as the stars rotate. 



Meibom et al 



fl2011bh published 



The above situation started to improve recently when 
rotation periods for stars in the 1 Gyr old cluster NGC 6811. Theirs is a comprehensive 
program that takes advantage of the existence of three open clusters within the field of 
view of NASA's Kepler mission, and thus it is expected that their work for NGC 6811 will 
be repeated for the cases of NGC 6819 and NGC 6791, which would provide calibrations of 
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gyrochronology relations at 2.5 and, possibly, 9 Gyr of age, respectively. 

In this paper we introduce an ongoing program devised to place constraints for 
gyrochronology relations in the regime of ages where the large majority of field stars fall, 
i.e., from one to several Gyr. We achieve this by studying samples of wide stellar binaries 
chosen so that one of their components is an evolved star (either turnoff, subgiant, or 
white dwarf), while their secondaries are regular MS stars of FGK types. Therefore, with 
the evolved primary potentially well suited to provide a reliable age for the systeno, a 
measurement of the rotation period of the MS secondary automatically provides a point 
useful for the calibration of gyrochronology relations at the corresponding age. 

The availability of large numbers of wide binaries with the above characteristics, both 
from existing as well as upcoming all-sky surveys and catalogs, suggests that these objects 
could in principle provide quite a large and varied pool of gyrochronology constraints. 
Indeed, being a representative sample of the mix of stellar populations that surround the 
solar neighborhood, wide binaries span a large range of stellar properties, including age, 
mass, and metallicity. For these reasons, a gyrochronology program based on these objects 
has the potential of populating the phase space of relevant variables (i.e., age, mass, and 
rotation period) more densely and continuously than it may be possible with only star 
cluster observations. 



^The assumption of a common age for the components of a wide bin ary should be a safe 



one. In the context of the recently proposed formation mechanism of 



Kouwenhoven et al. 



( 120101 ). where wide binaries are formed during the dissolution phase of young star clusters 
when pairs of initially unassociated stars get close in phase space and become bound, the 
maximum difference possible between the ages of the two components cannot be larger than 
the age of the dissolving young cluster (on the order of 10^ — 10^ yr), thus small, when not 
negligible, in comparison with the present Gyr-scale ages of these systems. 
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The first step in our program, therefore, involves tlie determination of stellar ages for a 
sample of wide binary components suitable of being dated. Today, precise stellar ages for 
individual stars can be readily computed for stars that either (a) began relatively recently 
to evolve away from the MS, such as turnoff stars and subgiants, or (b) did this a long time 
ago and are now in the white dwarf (WD) phase of their evolution. For recently evolved 
stars, reliable ages are routinely obtained via the use of adequately-chosen isochrones, 
provided stellar parameters such as effective temperature (Tes), metallicity ([Fe/H]), surface 
gravity {logg), and distance are known with precision. In the case of a WD, if it is of a DA 
type (i.e., with a hydrogen-dominated atmosphere), its cooling time (tcooi) and mass can 
be derived from its T^s and logg (measured from its spectrum) and the use of appropriate 
cooling models. Then, using empirically-calibrated initial-to-final mass relations one can 
derive the mass and hfeti me of the progeri i tor star, which, when added to tcooi provides 



the total age of the WD fIZhao et al.l 12011 



IGarces. Catalan, fc Ribas 



201l[ ). Thus, wide 



binary systems where one of the components is of either of these types is suitable to be 
age-dated, and this is what we plan to take advantage of in order to assign ages to their 
FGK components, which could then serve as gyro chronology constraints. 

Once their ages have been measured, the remaining step for obtaining gyro chronology 
constraints requires the measurement of the rotation periods of the MS components of these 
pairs. Rotation periods of individual stars are typically obtained by monitoring the small 



modulation o f the star's brig. 



rotatei 



(e.g. 



itness as spots come in and out of the line-of-sight as the star 



Hartman et al. 



20111 ). The monitoring can also be done spectroscopically, by 



^Rotation periods obtained in this way are thus independent of the inclination of the 
rotation axis, i.e., there is no sini ambiguity as in the case of the measurement of the 
rotation velocity from the broadening of the star's spectral lines. Transforming from a 
rotation velocity to a rotation period would suffer, additionally, from the uncertainty on the 



- 8 - 



following the intensity v ariation of the emiss i on at the cores of the Ca II H and K 
as the star rotates (e.g.. 



Hall. Lockwood. &: Skifl 



Vaughan et al. 



1981 



Baliunas et al 



1983 



Cincunegui et al. 



mes 



20071 : 



20071 ). We have not been able to find published rotation periods 



for any of the MS components of the pairs presented in this paper. These targets are bright 
stars {5 < V < 10 mag), and typical stellar rotation periods at such old ages are of the order 
of a month and longer. Two space missions, the MOST and CoRoT satellites, are very well 



posed for this kind of work, routinely achieving preci sions of a few mil 



smal 



2010 



er on stars similar to those of our program (e.g . 



Miller-Ricci et al. 



imagnitudes and 



2008 



Siwak et al. 



Strassmeier et al. 



2010 



Csizmadia et al. 



201ll ). From the ground, the only possible 



way to achieve this is using almost-dedicated telescopes with small apertures, and for this 
purpose we started to take advantage of instrumentation of that kind at the Observatorio 
Docente of Universidad Catolica in the outskirts of Santiago, Chile. This part of the 
program is currently in progress for an initial sample of the systems presented here, and 
will be reported upon completion in forthcoming papers. 

In this paper we describe our selection of an initial sample of wide binaries suitable 
for the purposes of this program and report on the ages of a number of these systems, 
obtained either from our own methods or directly from the literature. Section |2] describes 
the selection of targets. In § 3 we describe our procedure to determine accurate stellar 
parameters for a subset of target stars to be used in the process of age determination, which 
is reported in § 4. In § 5 we summarize our program and initial results. 



actual size of the star. 
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2. Sample Selection 

There exists a number of published wide binary catalogs that can be used to select 
systems suitable for our gyrochronology program. These catalogs have been assembled from 
a variety of parent surveys and following a variety of selection criteria regarding angular 
separation limits, proper-motion cuts, stellar type, and even distances. Since the crucial 
assumption lying at the core of this program is that a number of properties that can be 
well measured for one of the components of a wide binary can be safely assigned to the 
other componenjfl, we need to maximize the chances of selecting genuine, truly bound pairs 
of stars, avoiding any possible contamination by unassociated pairs as much as possible. 
Therefore, when searching for useful wide binaries for our program the best sources will be 
those that provide the largest amount of information (position, photometry, kinematics) 
that can be used to assess the evidence for binarity on any given pair of stars. Moreover, 
given that useful gyrochronology constraints require the determination of reliable ages, the 
best targets for our program will be those for which stellar parameters can be measured 
with high precision. 

The publicly available wide binary catalogs that best satisfy the above requirements 
are those of Chaname & Gould (2004; hereafter CG04), Gould & Chaname (2004; hereafter 
GC04) and Lepine & Bongiorno (2007; hereafter LB07). Although the literature offers 
a few catalogs with larger numbers of binaries than these three, most nevertheless suffer 
from shortcomings that would negatively impact a sample aimed for gyrochronology: 



large contamination 



jy the chance alignment of unassociated pairs (e.g., the catalog of 



Sesar. Ivezic. fc Juric 



2008 



based on SDSS data), faint sta rs comprising the large majority 



of objects in the catalog (e.g., the SLoWPoKES catalog of iDhital et al. 



2010. also based on 



^For close binaries, of course, this statement might not always be true, as their components 
may have interacted with each other and thus affected their normal evolution. 
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SDSS data and which, moreover, was aimed by construction at very low-mass, late-K and 
M stars), and even a significant fraction of systems in the catalog expected to be not bound 



anymore (IShaya fc Oiling 



20111 ) 



In contrast, the CG04, GC04, and LB07 catalogs only contain pairs with high 
probability of being genuinely bound systems, satisfying stringent requirements not only of 
kinematical nature (i.e., both stars displaying a common proper- motion), but also having 
luminosities and colors consistent with the two stars sharing the same age and chemical 
composition. Moreover, unlike those based on SDSS data, the CG04, GC04, and LB07 
catalogs were assembled from surveys of high proper-motion stars, where the probability 
of two stars being close in phase space and moving at such high velocities but not being 
associated is already very small. Finally, having been selected from a high proper-motion 
survey, the stars in these three catalogs are among the brightest on the sky, and hence 
the best suited to provide precise stellar parameters through high-resolution spectroscopic 
work. 

Before proceeding to the selection of targets, we briefly discuss t he impact on our 



program of undetected higher-order rn ultiplicity among wide binaries (jCorreia et al. 



2006 



Tokovinin. Hartung. fc Hayward 



20101 ). If the evolved star of a pair in our sample has an 
unresolved, undetected close companion, in principle the stellar parameters measured for 
the program star may be affected by light from the undetected neighbor, and thus the age 
derived for the system. However, as long as such hypothetical companion is signiflcantly 
fainter than the primary star, for high signal-to-noise observations such as those we deal 
with in §[3] the impact of an unresolved companion would be minimal. Moreover, since we 
are working in the visible part of the spectrum, the effect is even more negligible. The 
impact of a cool, close companion on the observed spectrum of the evolved star would be 
more important in the infrared and at longer wavelengths. On the other hand, if the star 



- 11 - 



with an unresolved companion is the FGK-type MS member of the pair, then there is a 
chance, depending on how close the unseen companion is, that the rotational properties 
of the program star have been affected by interactions with such companion, and thus its 
surface rotation rate may not reflect the processes behind gyrochronology. This is a worry 
in common with gyrochronology programs based on star clusters, and indeed many of them 
include a campaign of spec troscopic monitoring in order to identify close binaries (e.g.. 



Meibom et al. 



2009 



2011al ). In conclusion, pairs with MS components having signs of a 



close companion should be avoided and are thus excluded from our program. 



2.1. Wide binaries with turnoff and subgiant components 

Given that isochronal ages are very sensitive to uncertainties in distance, we restrict 
our selection to binary systems for which at least one of the components is an Hipparcos 
star. This is the case of all entries in the GC04 and LB07 catalogs, and all pairs in CG04 
with this characteristic are by construction already listed in GC04. Therefore, in order to 
select wide binaries with recently evolved components, we focus on the GC04 and LB07 
catalogs and restrict our search to pairs with better than 3cr parallaxes (vr/cr^ > 3). 

As a first step, for all entries in the GC04 and LB07 catalogs, we make the distinction 
between pairs in which both components are Hipparcos stars with independent parallax 
measurements, and wide binaries in which only one of the components is in Hipparcos. This 
is simply because having independent trigonometric distances to both stars in any given 
pair provides an important extra criterion for assessing the true binary nature of the stellar 
pair, thus placing these cases in a different category. 

We proceed to identify systems with evolved components based on the comparison 
with the well known Hipparcos CMD obtained with B — V colors and l^— band absolute 
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magnitudes. This is illustrated in Figures [T] and |2l corresponding to the two distinct 
categories we defined above according to the existence of independent Hipparcos 
measurements for both or just one of the binary components, respectively. The small 
background dots in the lower plots are Hipparcos stars within 100 pc of the Sun and 
with better than 5cr parallaxes, obtained from the Vizier Serviccl. They clearly show the 
population of field stars we are after, i.e., evolved stars leaving the MS on their way to 
become red giants. The top panels of these two figures show a CMD based on the V — J 
color, which is the reference color in the GC04 and LB07 catalogs. The small dots in the 
V — J CMDs are a subset of the LSPM-North p roper-motion catalog that includes stars 



20051) • Since not listed in GC04 



with good distances within 33 pc from the Sun (iLepind 
and LB07, optical B — V colors for all stars in the catalogs were also obtained from Vizier 
queries to the Hipparcos catalog. 

Wide binaries with recently evolved primaries were selected by defining the area 
contained within the dashed lines in the lower CMDs of Figures [1] and [21 This area was 
designed so that it encompasses the largest fraction possible of turnoff stars and subgiants, 
and at the same time attempting to avoid serious contamination from MS stars on the blue 
side. Given the smaller density of stars on the red side, we were less restrictive on that end 
and extended the selection box so that it includes the base of the red giant branch. 

From the pool of wide binaries in GC04 and LB07 with both components in Hipparcos, 
we selected all pairs for which any of the two stars fall within the search area defined above, 
finding 6 with MS components bluer than B — V ^ 1.40, approximately the boundary 
between K and M dw arfs. Since our source catalogs are based on the original Hipparcos 



reduction (ESA 



19971 ). we performed this initial search using parallax measurements 



from that database. Repeating the same exercise using the parallax measurements from 



"^http:/ /vizier. u-strasbg.fr/viz-bin/VizieR 
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the revised Hipparcos reduction ( Ivan LeeuwenI 120071 ) produced two more pairs, which we 



included in the sam ple. An additional pa ir, HIP 10305/HIP 10303, long known to be a wide 



binary system (e.g.. 



van der Bergh 



19581 ). was part of a sample on a project different from 
that described in this paper, but, finding it suitable for gyro chronology, it was included in 
our sample. The 9 pairs selected in this way are shown in Figure [T] and listed in Table [H 
The parallaxes reported in this Table and throughout this work are from the 2007 Hipparcos 
reduction. 

The second pool of wide binaries from GC04 and LB07 are those with only one of their 
components in Hipparcos. The selection of systems with a recently evolved component from 
this group is performed exactly as done for the first group above, only that this time the 
secondary is assigned the parallax measured for the primary. Since we could not find B — V 
measurements for all the secondaries in this group, we only required for these stars that 
V — J < 3.0, which approximately marks the boundary between K and M dwarfs. The 10 
targets found from this group are illustrated in Figure |2] and listed in Table |2l Note that 
one of these pairs is composed of a turnoff/subgiant primary and a WD secondary. While 
not useful for gyrochronology purposes, we decide to include it in our program because it 
may at some point serve us to compare the age of the system as derived from two different 
methods. 

As will be seen in §[3l the radial velocities and metallicities that we derive for an 
initial subsample of the targets selected in this section confirms that these are indeed 
binary systems. The mean [Fe/H] difference is 0.013 ± 0.084, i.e. consistent with zero, 
and the radial velocities agree within ~ 2 km/s (the mean difference for 10 pairs is 
— 1.05 ± 1.76 km/s). The exception are three pairs showing radial velocity differences of 
3-4 km/s (HIP 94076/HIP 94075, HIP 115126/NLTT 56465, and HIP 114855/NLTT 
56278). For truly bound wide binary sytems, such large differences between the radial 
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velo cities of their components can only be explained by detected orbital motion (see figure 



1 in 



Chaname fc Gould 



20041 ). or else if one of the components has a close, undetected 
companion that induces periodic radial velocity variations. Indeed, this is the case of HIP 
94076, a known spectroscopic binary, and we plan to obtain new epochs of spectroscopy for 



the other two cases. 



2.2. Wide binaries with white dwarf components 



The CG04 catalog contains 82 pairs with a WD component, and GC04 contains 20 
pairs composed of a WD and an Hipparcos star. These numbers, however, include all types 
of WDs as well as all types of MS components, while we are interested here only in pairs 
composed of a DA WD and a MS star of FGK type. 

A first selection and age-dating of pairs from CG 04 and GC04 containi i ig W' D 



Garces. Catalan, fc Ribaj (120111 ). who 



components of DA type was recently performed by 
took advantage of such systems in order to obtain the ages of a number of GKM stars 
aimed to calibrate the time evolution of high-energy emissions associated to stellar activity 
in low-mass stars. Although the majority of wide binaries in their sample contain M dwarf 
components, there are at least 7 systems with GK components that could be useful for 
the purposes of gyro chronology. We list these in Table [3] along with relevant photometric, 
astrometric, and kinematic data, as well as the ages derived by their work. 

Addi tionally, we obtain 11 more pairs composed of a WD plus a GK main sequence 



star from 



Zhao et al. 



(120111 ). who investigated the time evolution of chromospheric activity 
levels in solar-type dw arfs belonging to wide binaries selec t ed fro m the original Luyten 
proper motion survey (iLuyteniil979i) and from iGiclas et al.l (119711 ). Relevant data for these 



pairs are also listed in Table [3l when available. 
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We also inspec ted the 21 pairs with WD components in the SLoWPoKES catalog 



(jPhital et al. 



20101 ). However, all stars in this sample are faint (both primaries and 
secondaries in the range r = 16 — 20 mag) and, moreover, all the MS companions of these 
WDs turned out to be M dwarfs. 



3. Stellar Parameters for an Initial Sample 



As a first step in our gyrochronology program, we set out to determine the ages of 
a number of the wide binary systems selected in §[2l The sample of wide bina r ies wi th 
WD components of § 12.21 has already been studied by lGarces. Catalan, fc Ribad (l201ll ) and 



Zhao et al 



(I2OIII ). and thus we simply take their resulting ages and list them in Table [3] 
along with all the relevant data for these systems. In what follows, we concentrate on the 
samples of wide binaries with turnoff or subgiant components. 

Since the ages of recently evolved stars are derived with the use of theoretical isochrones, 
we first need to obtain sufficiently precise measurements of the stellar parameters (Tefj, 
[Fe/H], and logg) of the evolved component in these targets. In principle, this is only 
necessary for the component of the binary that is best positioned to provide information 
on the age of the systeno- Due to a number of reasons, however, we attempt to obtain 
isochronal ages for both components of our pairs, when possible. First, the procedure of 
age-dating via isochrones not only provides ages but can also be used to better characterize 
the stellar properties and evolutionary state of the stars being studied. This is important in 



^While evolved stars have typically been the best type of star for this, we note that today 
very high-precision stellar parameters can be derived for solar analog stars, and those can 



in turn be used to derive reasonably 



mam-sequence (e.g., 



Baumann et al. 



good isochrone ages, even for unevolved stars on the 



2OI0I ). 
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our context because the entire idea behind gyro chronology makes sense only for stars on the 
MS, whose rotation periods respond to their spin-down due to stellar activity. Therefore, 
we need to rule out as best as we can the possibility that some of the MS secondaries of the 
wide binaries selected in §[2] may have started themselves to evolve away from the MS, at 
which point the physical expansion of the star would affect its surface rotation rate and 
thus invalidate the system as a useful constraint for gyro chronology. Second, even for the 
cases where the MS components of the selected binaries are indeed FGK stars on the MS, 
we want to explore the possibility of obtaining a better constraint on the age of the system 
by forcing the procedure of isochrone fitting to consider the two stars as coeval. Third, for 
the cases where both stars are able to provide independent isochronal ages for the system, 
we would like to explore closely those pairs that produce inconsistent results, if any. Fourth, 
independent determinations of the metallicities of both components of our wide binaries 
automatically provide us with a check on our measurement errors, and may also serve as 
an additional test for the binary nature of our pairqj. Therefore, for the measurement of 
stellar parameters prior to age determination, we target both components of our selected 
wide binaries. 



Tn the context of the scenario of 



Kouwenhoven et al. 



(120101 ) for the formation of wide 



binaries, there is the possibility that the components of a genuine, gravitationally bound 
system may truly have slightly different metallicities, which would correspond to metallicity 
variations or gradients across the parent young star cluster. Additionally, it is becoming 
i ncreasingly more clea r that small abundanc e changes could be produced by planet formation 



( iMelendez et al. 



2009 



Ramirez et al. 



systems have been detected too (jPesidera et al. 



201 of), and small differences betwee n stars in binary 



2004 



Ramirez et al. 



20111 ). Thus we expect 



the metallicities of the two stars in a binary system to be similar but not necessarily identical 
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3.1. Observations and Data Reduction 



Most of our target stars accessible from 



with the MIKE spectrograph (iBernstein et al 



;he So uthern hemisphere were observed 



20031 ) on the 6.5 m Clay Telescope at Las 



Campanas Observatory on 2010 September 21-22 and 2011 January 4. We used a narrow 
slit (0.35"), which delivers data with spectral resolution R = A/AA 65 000 (at A ~ 6000 A) 
and the standard setup that allows complete wavelength coverage in the 3400-9100 A 
spectral window. These spectra were reduced using the CarnegiePython pipelinejzl which 
employs multiple bias and flat-field frames to correct for instrumental imperfections and 
ThAr lamp exposures taken throughout each night for wavelength calibration, in addition 
to co-adding multiple exposures of the same object. The signal-to-noise ratio (S/N) of our 
reduced spectra (per pixel) varies between about 100 and 600 at A = 6000 A with a median 
of about S/N=400 (at A = 4000 A the median S/N is about 150). 

Reflected sun-light spectra from the asteroid Hebe were acquired on 2010 September 
21 for solar reference. Observationally, asteroids behave like point sources, thus making 
their data acquisition and reduction identical as that for the rest of our targets. This is 
not the case of scattered sky-light or Moon observations, which are sometimes also used for 
solar reference. The use of asteroid observations allows a more precise differential analysis. 

The radial velocities of the stars were estimated from the Doppler shifts of spectral 
line cores of hundreds of Fel fe atures. We used the rest wavelengths measured in the 



laboratory by 



Nave et al. 



(Il994l ) and measured core wavelengths in the observed spectra 
by fitting a parabola to the 7 data points closest to the flux minimum of each line. The 
internal precision of our radial velocity measurements is about 0.35 km s^^. However, 
systematic errors due to core wavelength shifts produced by surface convection are of order 



^http://obs .carnegiescience.edu/Code/mike 
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0.5 km s ^ (e.g. 



Gray 



2009 



Ramirez et al 



20091 ). The use of cross-correlation with radial 



velocity templates would not remove entirely this error because of the wide range of stellar 
parameters of our sample and the fact that the impact of granulation is still poorly known 
even for standard stars. Nevertheless, the accuracy of our radial velocities (~ 0.5kms~^) 
is sufficient for our purposes. Our observing log and radial velocities derived are listed in 
Table H 



3.2. Spectroscopic Analysis 



The fundamental atmospheric parameters T^s, logg, and [Fe/H] of a star can be 
estimated using a variety of techniquesjf] In addition to employing only the observed 
spectra, photometric data as well as trigonometric parallaxes can be used to constrain 
one or more of these quantities. Here we describe the techniques used in our work. The 
atmospheric models ado pted are from the MARCS grid of standard chemical composition 



(IGustafsson et al. 



20081) o The curve-of-gr owth analysis was made using the 2010 version of 



the spectrum synthesis code MOOG (e.g.. 



Sneden 



19731)1^° 



We started by using the standard iron line spectroscopic approach, which forces 
excitation and ionization balance of iron lines. A first guess of the parameters is made 
and iron abundances are computed for a number of neutral (Fel) and singly ionized 
(Fell) iron lines. The parameters are then fine-tuned to remove any correlation between 
iron abundance and excitation potential (EP) of Fel lines (therefore forcing excitation 



^We use the standard notation: [Fe/H] = Ap^ — Ap^, where Ap^ = ^og{Npc/Nii) + 12 and 
Ax is the number density of X atoms in the stellar photosphere, 
^http:/ /marcs. astro. uu.se 

^°http:/ /www. as. utexas.edu/~chris/moog. html 
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balance) and to minimize the difference between the mean iron abundances inferred from 
Fel and Fell hnes separately (thus achieving ionization balance). Simultaneously, the 
correlation between Pel abundance and line strength is controlled (i.e., the correlation is 
minimized) with the microturbulent velocity parameter (vt). Fig. [3] shows an example of 
the e nd product of t his pr ocedure. The iron line-list adopted (including atomic data) is 



from 



Asplund et al. 



(120091 ). who made a careful selection of unblended features for their 
solar abundance analysis. The abundances used in this procedure are differential, on a 
line-by-line basis, using the solar abundances inferred from our solar (Hebe) spectrum as 
reference. We adopted vf = l.Okms"^, although the exact value of vf has a minor impact 
on our results. 

Errors in the derived parameters are estimated from the uncertainty in the abundance 
versus EP slope (for Tgfj) and the line-to-line scatter of the mean Fel and Fell abundances 
(for log (7). Since we force the EP slope to be zero, a slightly positive (negative) EP slope 
implies a Teg too low (high) by a certain amount. We use the ATes amount that corresponds 
to an EP slope of ±1(t, where a is the error of the zero slope when using the adopted Tgff. 
For the error in log (7, we consider the maximum and minimum logg values such that the 
mean Pel — Fell abundance difference is consistent with zero within the la line-to-line 
scatter as the upper and lower limits of the derived logg. 

A second set of parameters was obtai ned using colors to deriv e Tgff. The recent 



metallicity-dependent color calibrations by 



Casagrande et al.l (I2OIOI ) for the following color 



indices: B — V,b — y,V — J,V — H,V — Kg, and J — Kg, were used. Observed magnitudes 



and colors of our samp 



(IMermilliod et al. 



1997 



e stars were taken from the General Catalogue of Photometric Data 



11 



and the 2MASS and Hipparcos/Tycho catalogs. We made sure 



that the adopted photometry was not blended (i.e., we excluded mainly old measurements 



11 



http: / / obswww.unige.ch/gcpd/gcpd.html 
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for the unresolved systems) a nd avoided uncertain 2M ASS photometry for the brightest 



stars. We did not apply the 



Casagrande et al. 



(120101 ) formulas to giant stars because 



their work is restricted to dwarf and subgiant stars. Errors in the photometry and the 
color-to-color Teg scatter provides us with an estimate of the error in Tcs. Photometric 
errors are simply propagated into the color— Tcfj relations to obtain the error in T^s for a 
given color. Then these errors are used as weights when computing the final T^s value. The 
Teff error is obtained using the formula for the sample variance, again using the T^s errors 
for each color as weights. 

Given a photometric T^s, the surface gravity was determined using two methods. 
The first one is the same as in the iron line analysis, i.e., log (7 is fine-tuned so that the 
mean iron abundances from Fel and Fell lines agree (ionization balance). In the second 
case we determi ne log q using the s tars' trigonometric parallax from the new reduction of 



Hipparcos data (Ivan Leeuwen 



20071 ) and theoretical isochrones. This method is described 
in detail in §|U Note that in this latter case it is not guaranteed that the iron abundances 
inferred from Fe I and Fe ll lines are the same. Moreover, either if log g is inferred forcing 
ionization balance or using isochrones, the Fe I abundances will in general show a correlation 
with EP. Thus, the line-to-line scatter of the iron abundances inferred using photometric 
temperatures will be larger than that obtained by forcing excitation and ionization balance 
of iron lines. However, this does not imply a superiority of one method over another; it 
simply reflects the nature of the different approaches to measure the stellar parameters. 

Finally, a third estimate of effective temperature can be obtained by analyzing the 
wings of Balmer lines, in particular Ha. The depth of these wings is highly sensitive to 
Teff and if the other stellar parameters can be constrained independently, very precise Tes 



values can be inferred from a minimization of observation rninus theoretica 



Ha line profiles (e.g. 



Barklem et al. 



2002 



Ramirez et al. 



2006 



Ramirez et al 



models of 



2OI1I). A 
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proper continuum normalization is required for this method to provide accurate effective 
temperatures. Our Ha hne profiles were normalized taking advantage of the smooth 
variation of the blaze function across spectral orders. Polynomial fits were used to trace 
the upper envelopes of spectral orders above and below the order containing the Ha line. 
They were then interpolated to trace th e continuum of the H a order. We used a grid of 



theoretical Ha line profiles computed by 



iBarklem et al. 



(]2002h . which is based o n MA RCS 



Barklem et al. 



((20001)0 A 



atmospheric models and the self-broadening theory developed in 

minimization routine allowed us to find the best model fits to our data and therefore 
to estimate T^e and its associated error. Fig. H] illustrates this technique. We obtained 



Teff = 5741 ± 4 K for our solar (He 



Teff inferred by; 



spectrum by 



Barklem et al. 



Kurucz et al. 



De) spectrum, in very good agreement with the solar 



([2OO2I ) using the Ha line from the very high quality solar 



19841 ). The fact that the Ha temperature of the Sun does 



not perfectly agree with the nominal value of 5777 K suggests that there is still room for 
improvement in the modeling of Balmer lines. We increased the Tefj values derived from our 
Ha analysis by 36 K as a first order correction, given that this brings the solar Teg up to its 
expected value. Similar to the case of photometric Tgfj, once the effective temperature was 
determined from the Ha line analysis, we derived \ogg values using ionization balance as 
well as isochrones. 

Note that all the methods described above require a previous knowledge of the 
atmospheric parameters, which are the quantities we want to derive. Thus, iterative 
procedures are necessary to obtain a final, self-consistent solution for the Teg, log (7, [Fe/H] 
set. We repeated all calculations as many times as necessary to guarantee that no more 
iterations are required to improve these solutions. 



12 



This grid of theoretical Ha line profiles is available online at 



http: / /www.astro. uu.se/ ~barklem 
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The Teff, log (7, [Fe/H] values derived using the various techniques described in this 
Section are listed in Table O 



3.3. Adopted Parameters 



All methods of atmospheric parameter determination have limitations. For example, 
the excitation/ionization balance of iron lines assumes local thermody namic equilibri um 



Asplund 



20051). 



(LTE), which is a useful but not necessarily correct assumption (e.g., 
Photometric temperature scales have uncertain zero points, and even though they seem to 
agree well with Tgfj scales based on interferometric observations of stellar angular diameter, 
the latter are typically based on the analysis of giant stars, whose surfaces are far from 
being static and well-defined, as it is assuined in the measurement of angular diameters 



[e.g., 



iKoesterke et al. 



2008 



Chiavassa et al. 



20101). Also, the modeling of Ha line profiles 



relies heavily on the theory of self-broadening, and different prescr i ptions may lead to 



Barklem et al. 



(I2OO2I ). Surface gravity 



significantly different results, as explored in detail by 
determinations based on ionization balance rely also on the assumption of LTE while the 
isochrone approach depends on st ellar structure and evolution calculations and ra ay be 



affected by statistical biases (e.g.. 



Pont fc Ever 



2004 



J0rgensen fc Lindegren 



2OO5I). 



It is therefore not surprising to find differences in the stellar parameters derived with 
different methods for the same object in Table O Our sample size is too small and it 
covers too wide of a region of stellar parameter space for us to clearly uncover systematic 
differences between the different methods. In any case, discrepant values typically point 
to systematic errors affecting differently each technique whereas values in good agreement 
suggest that the impact of these errors is relatively small. Thus, by taking the average value 
of all the measurements reported in Table [5] for each object we improve upon values that 
are only slightly affected by systematic errors while assigning realistic error bars to stars for 
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which their derived parameters hkely suffer from severe systematic errors. 

Our adopted parameters are obtained as the weighted average of all measurements 
available to each star. However, effective temperature errors lower than 35 K were set to 
35 K before averaging. Similarly, the lowest adopted logg and [Fe/H] errors were 0.03 and 
0.04 dex, respectively. This was done to avoid giving too much weight to a particular value. 
Small internal errors can in many cases be the result of numerical artifacts that do not 
reflect the true errors of the measurement. By adopting reasonable minimum errors we 
obtain more reliable, and therefore accurate, average values. 

The errors adopted are the linear sum of the standard error and the sample variance. 
The first term takes into account the random errors of each measurement assuming that 
there are no systematic differences between the various values. The sample variance, on the 
other hand, is sensitive to systematic differences. Thus, by adding the sample variance to 
the standard error a more realistic error bar is obtained. 

Our adopted parameters and their associated errors, computed as described in this 
Section, are given in Table [61 



4. Stellar ages 

We compute the ages of our initial sample of recently evolved wide binaries by making 
use of theoretical isochrones. In this technique, the star under study is placed on a CMD 
and its location compared to theoretical predictions of stellar evolution. Isochrone points 
:o the observed stellar parameters are then used to derive the mass and age of the star 



close 
(e.g., 



Lachaume et al. 



19991 ). While the location of any star on the CMD is determined just 
by its absolute luminosity My and surface temperature Tcs (or color), many combinations 
of stellar parameters (}ogg, [Fe/H], mass, and age) can go through that point, so that in 
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order to infer the age of the star some of these parameters must also be given. Typically, the 
fundamental atmospheric parameters T^s, logg, and [Fe/H] are measured from the star's 
spectrum, photometry, or a combination of both (§[3]), and they are enough to determine 
both the age and mass of the star via isochrones. In our case, given that for some stars 
different methods to estimate logg produced discrepant answers (§ l3.3p . we avoid this 
parameter and instead use the luminosity of the star. My (i.e., determined by its apparent 
magnitude and parallax). 

In order to determine which isochrone fits best any particular star, we build the 
probability distribution of the stellar age by computing the likelihood that any given 
isochrone passes near the corresponding set of stellar parameters, accounting for the 



uncertainties on those parameters. The age probability distribu 
"G function", and the procedure is widely used nowadays (e.g., 


tion is also callec 


. the 


Lac 


haume et al. 


1999; 


Reddv et al. 


2003; 


Nordstrom et al. 


2004; 


AUende Prieto et al. 


2004; 


Pont & Ever 


2004) 



Assuming the errors in our stellar parameters {(JT^g, <y[Fc/H], and <Jm^i) have Gaussian 
probability distributions, the likelihood for a point in a given isochrone can be written as 



P(Teff, [Fe/H], Mv) ocexp 



r-(ATeff)21 




exp 


L J-eS -> 



-(A[Fe/H])^ 

[Fc/H] 





'-(AMv)2" 


exp 


[ J 



where ATgfj, A [Fe/H], and AMy are the differences between the measured stellar parameters 
and those corresponding to the point in the isochrone under consideration. The integral 
of this likelihood over all the parameter space of the set of isochrones then gives the age 
probability distribution 



P(age) = j P{T,s, [Fe/H], My) dT,s d[Fe/H] dMy. (2) 
In practice, we only integrate over the volume defined by three times the la uncertainties 
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from the measured stellar parameters, which we verified already accounts for most of the 
contribution to the probability distribution from the entire set of isochrones. Being the 
most likely value, we then adopt the peak of the G-function as the age of the star. 

The choice on how to determine the exact location of the peak of the G-function has a 
small effect on the adopted age, especially given the size of our la errors, typically larger 
than ~ 0.5 Gyr. Therefore, we do not implement at this point anything more sophisticated 
than just adopting as our age the center of the bin where the peak occurs, and leave for a 
later stage in this program any refinement, if at all needed. Due to the discrete nature of 
the procedure of isochrone fitting, this lack of "smoothing" of the G-functions leads to the 
appearance of spuriou s structure in the distribution of ages of large samples of stars (e.g.. 



Nordstrom et al. 



20041 ). so a more sophisticated approach to determine the exact location 



of the peak is relevant for population studies, which is not the case here. 

As an illustration. Figure [5] shows the resulting age distributions and la uncertainties 
around the most likely age for the pair HIP 115126/NLTT 56465 (solid lines). We normalize 
the distributions so that the area below is equal to 1. The adopted errors are computed 
from the cumulative function of the age probability distribution (dashed lines), assuming 
the latter is well approximated by a Gaussian. However, as is clear from the bottom panel 
in Figure El this G-function is not always close to a Gaussian, and thus we adopt different 
la errors to both sides of the peak. Our la lower and upper limits mark the age interval of 
cumulative probability between 16% and 84%. 



Using this procedure we compute the ages of all stars (both primaries and secondaries 



2001 



Demarque et al. 



2004h 



in our initial sample. We use Yale-Yonsei isochrones ( lYi et al. 
sampled with constant steps of 0.1 Gyr in age and 0.02 dex in [Fe/H], covering ages from 0.1 
to 15 Gyr. In constructing the age probability distribution, we experimented with different 
bin sizes and found that 0.5 Gyr/bin was an optimal choice and smaller bins did not change 
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the resulting ages. For the case of binaries with both components in Hipparcos, we adopt 
for both stars the primary's parallax, which for all pairs in this category is the one with the 
better measurement (i.e., smaller error). The resulting ages and uncertainties for all stars 
that could be observed with enough signal-to-noise are listed in Table [7] 

Following the same steps leading to the derivation of the isochronal age via equations [T] 
and [21 one can derive an isochronal log g through the construction of the gravity probability 
distribution. Our results for this isochronal logg are listed in Table |6] (labeled as log^fiso), 
where they can be compared to the average of the resulting logg derived from all other 
methods. Since the logf^iso are obtained from our final Tgfj and [Fe/H] for each star, as 
listed also in Table [6l they are our preferred values, above the log g obtained as the average 
of different methods in Section [31 

Next we consider the interesting possibility of improving the ages of those systems in 
our sample for which the two components were able to provide an independent age. In these 
cases, since both components contain some information on the age of the system, instead 
of treating them independently we can treat them simultaneously by forcing the procedure 
of isochrone fitting to consider the two stars as coeval. This can be done by computing the 
age probability distribution of the binary system, which would be given by 



where Pa and Pb are the likelihoods for components A and B, respectively, and given 
by equation [H While equation [31 gives the same weight to both stars, one could assign 
different weights based on a number of criteria, but we choose to not attempt this. Figure M 
illustrates the case of HIP 115126/NLTT 56465, the same pair shown in Figure [5l The net 
effect of considering the two components simultaneously was to shift the entire G-function 
of the evolved primary towards younger ages, closer to the broad peak of the secondary's 




(3) 
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G-function. The nominal precision of the resulting combined age (i.e., the size of the la 
region around the peak), however, remains esentially the same as that for the evolved 
primary alone. 

All the ages derived up to this point come directly from the adopted set of isochrones 
and the fitting procedure detailed above. The use of a different set of isochrones will, in 
principle, change our results. Moreover, isochrone ages are known to suffer from a number 
of systematic biases arising, for example, from the fact that, due to the different timescales 
characteristic of different evolutionary stages, isochrones do not evenly map the relevant 
parameter space of stellar parameters. There is a signif icant literature on the subject of such 



age 



) iases and the available ways t o correct for them (jPont fc Eyer 



2004 



2004 



Nordstrom et al. 



J0rgensen Sz Lindegren 



20051 ). These systematics are important when dealing with 



large samples of stars and attempting to extract conclusions of a statistical nature, 
such as those related t o the age-metallicity relation in the solar neighborhood (e.g.. 



Casagrande et al. 



201ll ). Since the main purpose of the present paper is to introduce our 



program for gyro chronology with wide binaries and demonstrate that precise ages (i.e., with 
reasonably small internal errors) can be achieved for them, any refining of the ages derived 
in this section will be attempted later, if at all necessary, and thus we do not implement 
any corrections for statistical biases or derive ages from other isochrone grids at this point. 

Nevertheless, in order to provide readers with a quantitative idea of the differences 
expected when accounting for statistical biases and isochrone grid choice, we compute the 
ages obtained by making use of the "web interface for the Bayesian estimation of stellar 
parameters" , mainta ined at the Osservatorio Astronomico di Padova^^l and described in 



da Silva et al. 



(120061 ). This publicly avai lable tool uses a dif ferent set of isochrones than the 



one used in this paper (namely, that by iGirardi et al. 



20001) and accounts for some of the 



13 



http: / / stev.oapd.inaf.it / cgi-bin/param 
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biases discussed above. The tool requires the selec tion of pr i ors for the Bayesian analysis 



Chabrierl (j200l[ ). and a constant star 



for which we choose an initial mass fonction from 
formation rate in the interval from 1 to 12 Gyr. The resulting ages and their uncertainties 
are listed in Table [7] next to our own determinations. 

In Figure [7^ we show the comparison between our isochronal ages (for the components 
of the wide binaries treated indep e ndent ly) and those obtained accounting for statistical 



biases following the 



da Silva et al. 



( 120061 ) prescription and a different set of isochrones. In 



most cases the ages are remarkably similar, even in the adopted un certainties. The mea n 



da Silva et al 



(120061 ) 



difference of the maximum likelihood ages, for the cases where the 
ages are at least three times larger than their la errors, is only 0.1 ± 0.6 Gyr, thus giving 
reliability to our isochronal ages and suggesting that inclusion of the statistical biases is not 
crucial for our purposes. 

We also tested our age determination method against the Bayesian implementation by 
Casagrande et al. (2011, hereafter Gil). They derived improved stellar param eters and ages 



Nordstrom et al. 



20041 ). We selected 



of all stars in the Geneva- Gopenhagen Survey (GGS, 
a small sub-sample of 47 GGS stars from Gil. These stars have ages with errors smaller 
than 16% if young (< 2.8 Gyr) and smaller than 20% otherwise (> 2.8 Gyr), as given by 
Gil. We used the effective temperatures and metallicities given by Gil along with the 
stars' visual magnitudes and Hipparcos parallaxes to derive ages using our isochrone-fitting 
program, i.e., we used the exact same input data employed by Gil. Figure [7)d shows the 
comparison of Gil and our ages for this small sample. The agreement between the two sets 
of ages is excellent, with a mean difference of only 0.03 ± 0.45 Gyr. A closer inspection of 
this Figure suggests very small systematic differences. For stars older than 7 Gyr, our ages 
are slightly younger (—0.48 ± 0.33 Gyr) whereas for stars younger than 4 Gyr they appear 
slightly old (although the mean difference is still consistent with zero: 0.2 ± 0.4 Gyr). In 



- 29 - 



this test we employed the ages from Cll derived using P adova isochrones. A sim i 
using the ages they derived with BASTI isochrones (e.g. 



Pietrinferni et al. 



2004 



ar test 



20061) 



gives almost the same result. Thus, we show again that statistical biases and the choice of 
isochrone grid appears to have a minor impact on our project. 



5. Summary and Conclusions 

We have started a program that takes advantage of wide binaries with evolved 
components and solar-type main-sequence companions in order to obtain constraints for 
the calibration of age-rotation relations (gyrochronology) at ages of several Gyr. We have 
mined published catalogs of wide binaries and assembled a sample of 38 wide binary 
systems with either turnoff/subgiant (20) or white dwarf (18) components that are well 
positioned to provide ages for these systems. For an initial subsample comprised of 15 of 
the binaries with turnoff or subgiant components we measured precise stellar parameters 
from high resolution optical spectra, and used these results to derive ages for the binary 
systems using theoretical isochrones. As judged from a one-to-one comparison with results 
obtained by state-of-the-art analysis that takes into account a number of known systematics 
of the isochrone technique, our own isochronal ages and the adopted uncertainties appear 
to be robust. We obtained the ages of 16 of our 18 systems with white dwarf components 
directly from the literature, having been derived from models of white dwarf cooling and 
empirically-calibrated initial-to-final-mass relationships. 

The ages of the 31 wide binary systems provided in this work range from 1 to 9 Gyr. 
At least 15 of these systems have la uncertainties of about 0.5 Gyr or less. Those for 
which the rotation period of the main sequence companion is suitable to be measured will 
thus provide useful constraints for gyrochronology relations, in a regime of ages where the 
only constraint currently available is that provided by the Sun. Observations aimed at the 
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measurement of the rotation periods of these stars have started and will be reported in a 
future paper. 

As larger and more varied samples of wide binaries start to become available thanks 
to ongoing and planned all-sky surveys, the potential of programs like the one presented 
in this paper for the purpo ses of gyro ch ronology will only increase. Taking advantage of 



the SUPERBLINK survey flLepine 



20081 ). we have already more than doubled the sample 



reported in Tables [T] to [3l and observations aimed to derive their ages are scheduled. With 
the help of increasingly improving photometric distances and stellar parameters, it may 
be possible today or in the near future to extend this work to the fainter database of the 
Sloan Digital Sky Survey. Finally, the unprecedented scale, depth, and precision levels of 
the data expected from upcoming projects like Gaia and LSST will provide large quantities 
of wide binar y systems suita. ble not only for gyrochronology but for a large range of other 



applications ( IChaname 



20071 ) 
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Fig. 1. — Selection of wide binary pairs with turnoff/subgiant components, both members 
having independent Hipparcos parallaxes (primaries as filled symbols, secondaries as open 
symbols). The 9 pairs in this category are shown here on two different color- magnitude 
diagrams and listed in Table [H The small background dots in the B — V CMD are Hipparcos 
single stars within 100 pc of the Sun and with better than 5a parallaxes, and in the V — J 
CMD they are single Hipparcos stars in the LSPM catalog and within 33 pc of the Sun 
(Lepine 2005). Primaries within the dashed area in the B — V CMD and with main-sequence 
secondaries bluer than B — V = 1.40 are selected for further inspection. 
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Fig. 2. — Same as Figure [H but for pairs with only one Hipparcos component. The 10 
pairs shown here are hsted in Table [21 Primaries within the dashed area in the B — V 
color-magnitude diagram and with main-sequence secondaries bluer than V — J ^ 3.0, 
approximately the boundary between K and M dwarfs, are selected for further inspection. 
The magnitudes of a few pairs have been randomized by a small amount in order to avoid 
clutter. Due to lack of colors in some cases, not all of the MS secondaries can be shown. Note 
in the upper panel the pair composed of a turnofF/subgiant primary and a WD secondary 
(HIP 18824/NLTT 12412, with the WD located at My ~ 14 and 1/ - J ~ 1). Even though 

+Hp>rpi is nn ThT^TC TV/TJ-l star in this nnir wp lfpip>r> if in niir r>rncrrnnn in naaf +Vip> is nf DA 



- 39 - 



hip9148: 5650 4.02 -0.144 0.72 





0.0 
-0.1 



-0.2 E- 

-0.3 i- 



|.^a.^o.-<^....... 



5000 5500 6000 6500 7000 

Wavelength (A) 



Fig. 3. — Measured iron abundance as a function of excitation potential (top panel), reduced 
equivalent width (middle panel), and wavelength (bottom panel) for HIP 9148. Open circles 
(filled squares) correspond to Fe I (Fe ll) lines. Dotted lines are linear fits to the data. The 
legend on the top panel shows the star name followed by the derived parameters: T^s, logg, 
[Fe/H], and Vt. The solid line is at solar metallicity. 
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Fig. 4. — Top panel: Observed (solid line) and theoretical (dashed and dotted) Ha line 
profiles for HIP 9148. The dashed line corresponds to the best fit T^s Ha model while dotted 
lines show profiles for T^s ± 200 K. Bottom panel: reduced of observation minus model 
differences as a function of model Tes (open circles). Only the regions not affected by narrow 
spectral lines are used in this computation. The solid line is a parabolic fit to the 7 points 
closest to the minimum, which allows a better minimization. 
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Fig. 5. — Age probability distributions (solid lines) for the components of a representative 
case (not the best, not the worst) of our initial sample of wide binaries with recently evolved 
components. The vertical dotted lines indicate the location of the adopted la errors to 
both sides of the peak of the distributions, obtained from the cumulative function of the age 
distribution (dashed lines). The evolved primary provides a significantly better constraint 
on the age than the MS secondary. 
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Fig. 6. — Age probability distribution assuming coevality for the pair shown in Figure El i.e., 
considering the two components simultaneously (see text in §|1]). Upper panel: independent, 
single-star age probability distributions for HIP 115126 (A, primary) and NLTT 56465 (B, 
secondary). Lower panel: Joint probability distribution (solid line), obtained by forcing the 
isochrone fitting to pass through the two components of the binary using equation |3l 



-43 - 



CD 
O 
if) 
O 



cn 
< 



14 
12 
10 
8 
6 
4 

2 


14 

12 

10 



CJ 

I 



o 



< 



2 - 

hi. 



-T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 

(a) Gyrochronology sample 



9 



+ 



+ 



+ 



+ 



(b) Small sample of GCS stars 
with very precise ages 




4 6 8 10 
Age (Gyr) - TW 



12 



14 



Fie. 7.- 
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and makes use of a d ifferent set of isochrones. The filled circles are stars considered robust by 

1 

da Silva et al.l ((20061) (age/dage > 3). The agreement is excellent, with only 3 — 4 significantly 



discrepant cases out of a total of 24 stars, (b) Comparison of ages from t he GCS survey for 
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(j201ll ) and our ages, obtained using the exact same input data from Cll. 



Table 1. Wide binaries with recently evolved primaries, both stars in Hipparcos. 



HIP 


(mas/yr) 


Mdbc 
(mas/yr) 


(mas) 


(mas) 


V 
(mag) 


B-V 
(mag) 


V-J 
(mag) 


Sp.Typei 


A0 
(arcsec) 


(mas/yr) 


NLTT/LSPM 


Comments 




9243 
9247 


137.5 

138.0 


125.9 
129.4 


7.81 
7.82 


0.79 
0.94 


8.27 
9.11 


0.88 
0.54 


1.64 
1.02 


KOIV 

F8 


36.0 


3.5 


6664 
6665 






10305 
10303 


378.2 
383.5 


-70.6 
-72.2 


21.71 
20.89 


1.67 
7.39 


5.65 
7.56 


0.55 
0.72 


0.60 
1.03 


F8V 
G4 


16.7 


5.5 


7323 
7321 






58240 
58241 


-174.0 
-180.9 


-5.3 
-5.4 


21.13 
8.35 


6.41 
6.35 


7.80 
7.83 


0.67 
0.64 


1.31 
1.23 


G3V 
G4V 


18.8 


6.9 


29039 
29041 


pre-MS star^ 




74432 
74434 


-594.9 
-594.9 


296.3 
283.5 


36.71 
21.49 


1.91 
4.98 


6.69 
7.71 


0.68 
0.74 


1.23 
1.47 


G5V 
G7V 


23.4 


12.8 


39601 
39603 


secondary within box 


1 


81991 
81988 


-215.7 
-219.3 


-257.7 
-255.0 


21.80 
23.06 


0.91 
1.73 


6.55 
10.31 


0.89 
1.04 


1.61 
2.01 


G5 
K3 


163.0 


4.5 


43483 
43481 




hl^ 
hl^ 

1 


94076 
94075 


51.7 
69.1 


186.1 
196.6 


19.11 
23.48 


1.57 
4.95 


6.70 
7.97 


0.64 
0.75 


1.15 
1.22 


GIV 
G5 


15.9 


20.3 


47474 
47473 


spectroscopic binary^ 
substellar companion 




101916 

101932 


322.2 
317.8 


21.2 

20.1 


33.20 

34.90 


0.82 
1.14 


5.07 

8.52 


0.70 

0.91 


0.96 

1.58 


G5IV 

K2V 


214.0 


4.5 


49632 

49639 






101082 
101166 


68 
68 


221 
220 


15.69 
14.51 


0.50 
0.70 


5.97 
8.68 


0.94 
0.64 


1.47 
1.13 


KOIII 
G5 


214.2 


1.0 


J2029+8105 
J2030-I-8108 






102532 
102531 


-30 
-1 


-198 
-204 


25.82 
26.35 


1.20 
2.40 


4.25 
4.97 


1.04 
0.49 


0.97 
2.91 


KlIV 
F7V 


9.1 


29.6 


J2046-I-1607E 
J2046-H607W 







Note. — 



(1) Obtained from the SIMBAD database. 



Table 2. Wide binaries with recently evolved primaries, only primary star in Hipparcos. 



HIP 


NLTT/LSPM 




MDEC 






V 


B-V 


V-J 


Sp.Typei 


Ae 


A/j, 


Gomments 






(mas/yr) 


(mas/yr) 


(mas) 


(mas) 


(mag) 


(mag) 


(mag) 




(arcsec) 


(mas/yr) 




9148 


6581 


214 


50 


11.15 


1.11 


8.26 


0.70 


1.26 


G3V 


41.5 


5 






6583 


219 


52 






13.07 


9.99 


2.29 










18824 


12415 


373 


-12 


19.35 


0.63 


6.80 


0.63 


1.17 


GIV 


64.1 


4 






12412 


370 


-12 






17.82 


9.99 


0.99 








white dwarf 


23926 


14574 


-73 


289 


18.69 


0.49 


6.75 


0.66 


1.24 


G3V 


10.0 


16 






14573 


-64 


276 


18.46 


0.65 


10.11 


0.67 


1.62 


GO 






HIP 23923 


32935 


17118 


123 


156 


8.19 


0.57 


8.12 


0.91 


1.74 


KOIV 


46.1 


2 






17117 


123 


154 






10.11 


0.55 


1.09 


GO 








114855 


56282 


370 


-17 


21.77 


0.29 


4.24 


1.11 


1.98 


KOIII 


49.6 


14 


giant; substellar companio^ 




56278 


370 


-17 






9.73 


1.06 


2.42 


K3 






1 


115126 


56466 


266 


-45 


47.35 


2.47 


5.20 


0.79 


1.38 


G8.5IV 


12.2 


32 






56465 


266 


-45 






8.88 


0.88 


2.84 


K2V 








18448 


J0356+6950 


190 


-168 


6.44 


1.20 


9.34 


0.93 


1.68 


KO 


24.7 


1 






J0356+6951 


189 


-167 






14.72 


9.99 


2.52 










6914 


J0129+7412W 


194 


-121 


17.89 


0.60 


7.26 


0.67 


1.23 


G5 


27.0 


42 






J0129+7412E 


183 


-80 






12.61 


1.17 


2.89 










20197 


J0419+1416 


82 


-216 


11.89 


1.10 


7.54 


0.99 


1.81 


G7in 


215.0 


14 






J0419+1419 


87 


-203 






13.12 


0.75 


2.74 








possibly M dwarf 



107759 J2149-t-1402S 159 -41 5.78 1.50 9.65 0.74 1.41 G5 13.3 9 

J2149-I-1402N 154 -33 99.99 9.99 9.99 



Note. — (1) Obtained from the SIMBAD database. 



I 

O) 

I 



Table 3. Wide binaries with white dwarf components. 



NLTT 


(mas/yr) 


MDEC 

(mas/yr) 


(mas) 


(mas) 


V 
(mag) 


B-V 
(mag) 


V -J 
(mag) 


Sp.Typei 


A6I 
(aJTcsec) 


(mas/yr) 


Age 
(Gyr) 


Comments 


1762 
1759 


207 
222 


-53 
-44 


9.52 


1.63 


16.59 
10.28 


0.30 
0.78 


-0.078 
1.475 


DA 
Kl 


28.8 


18 


, ,-1-0.15 
-■-^-0.14 


HIP 2600 


29967 
29948 


-299 
-292 


-328 
-313 


22.94 


1.63 


17.26 
9.96 


9.99 
0.99 


1.19 
1.89 


DA 
K4 


202.8 


15 




HIP 59519 


13599 
13601 


228 
233 


-155 
-149 


56.02 


1.21 


15.94 
8.42 


0.65 
1.10 


1.34 
2.47 


DA/DC 
K7 


123.9 


8 


5 4+2-® 


HIP 21482 


55288 
55287 


548 
551 


-57 

-49 


28.72 


1.30 


16.50 
8.03 


0.42 
0.66 


0.87 
1.22 


DA 
G7 


41.8 


9 




HIP 113231 


44348 
44344 


-92 
-88 


252 
259 


18.18 


1.35 


17.50 
11.46 


0.10 
1.02 


9.99 
2.34 


DA/DF 
K7 


28.7 


8 




phot pix 


7890 
7887 


-116 
-119 


-145 
-155 


22.22 


2.50 


17.39 
9.84 


-1.19 
1.01 


1.24 
1.90 


DA 
K3 


40.5 


10 


5.0t?-i 


phot pix 


1374 

1370 


72 
74 


-205 
-206 


8.00 


0.65 


16.22 

12.90 


0.36 
1.70 


0.17 
2.16 


DA 

K6 


59.2 


2 


9 Q-H.6 


phot pix 


LP378-537 
BD+23 2539 


-101 


79 


56.7 




16.2 
9.7 


0.7422 


1.252 


DA 
KO 


20 




^•"-0.1 




LP786-6 
BD-18 2482 


-190 


33 


29.49 




15.7 
12.8 


-0.118 
1.0362 


-0.552 


DB 
K3 


31 




■'•>J_o.6 


NLTT 20260 


40 Eri B 
40 Eri A 


-2228 
-2240 


-3377 
-3420 


200.62 


0.23 


9.5 
4.4 


0.11 
0.82 


-0.349 
1.397 


DA 
Kl 


82 


45 


5 0+1-1 

'^•"-i.o 


NLTT 12868 
NLTT 12863 


L481-60 
CD-37 10500 


-415 
-415 


-214 
-215 


65.6 
65.13 


0.40 


12.8 
6.0 


0.30 
0.72 


1.051 


DA 
G7IV 


15 


1 


^•-■■-0.2 


NLTT 41169 
NLTT 41167 


L182-61 


-46 


-316 






14.1 


-0.09 


-0.214 


DB 




1 


1 4+1-0 


NLTT 16355 



Table 3 — Continued 



NLTT ^RA /^DEC T ^Tv V B - V V - J Sp.Type^ AO Afi Age Comments 

(mas/yr) (mas/yr) (mas) (mas) (mag) (mag) (mag) (arcsec) (mas/yr) (Gyr) 



CD-59 1275 


-46 


-316 


26.72 


0.29 


6.4 


0.59 


1.024 


GO 


41 






NLTT 16354 


CD-38 10980 
CD-38 10983 


76 
72 


1 

5 


76.00 
78.26 


2.56 
0.37 


11.0 
5.4 


-0.14 
0.597 


-0.548 
0.966 


DA 
G5 


345 


6 




HIP 80300 
HIP 80337 


LHS300B 
LHS300A 






32.3 




17.8 
13.2 






DB 
K 










LP592-80 
BD-1 469 


253 


-60 


14.89 


0.84 


17.2 
5.4 


1.04 


1.966 


DA 
KlIV 


49 






HIP 15383 


G216-B14B 
G216-B14A 










15.5 
12.0 




2.864 


DA 


23 




-'-•^-0.2 


B-band mag 


G116-16 
BD+44 1847 


38 
34 


-296 
-279 


34.6 
19.0 


4.00 
1.04 


15.4 
9.0 


0.21 
0.68 


0.355 
1.315 


DA 
GO 


1020 


17 


1 6+°-2 
'-■"-0.2 


NLTT 21338 
NLTT 21283 


G273-B1B 
G273-B1A 


169 


7 






16.4 
11.0 






DA 
G 


36 




o 4+0.8 


B-band mag 



Note. — (1) Spectral types are taken preferentially from the original source, i.e., either iGarces, Catalan, fc Ribaj 1 2011 ) or 
SIMBAD database. 



Zhao et al 



I 2011 ). or else from the 



Table 4: Observing Log and Derived Radial Velocities 



Object 


UT date 


UT time 


Exp. time 


S/N 


Radial velocity 




(yyyy-mm-dd) 


(hh:mm:ss) 


(seconds) 


A = 6000 A 


(kms~^) 


Hebe 


2010-09-21 


04:10:27 


1004 


289 




HIP9148 


2010-09-22 


03:38:02 


3648 


420 


17.13 ±0.25 


HIP9243 


2010-09-21 


04:47:48 


2670 


447 


7.84 ± 0.24 


HIP9247 


2010-09-21 


05:38:55 


3250 


314 


7.26 ± 0.50 


HIP10305 


2011-01-04 


01:43:37 


138 


479 


-3.81 ±0.78 


HIP10303 


2011-01-04 


01:05:54 


843 


505 


-4.90 ± 0.25 


HIP18824 


2010-09-21 


08:48:35 


2081 


111 


10.64 ±0.31 


HIP20197 


2010-09-22 


09:17:40 


405 


437 


-18.58 ±0.31 


HIP23926 


2010-09-21 


06:48:54 


435 


409 


44.30 ± 0.30 


HIP23923 


2010-09-21 


07:08:18 


4606 


180 


45.27 ±0.22 


HIP32935 


2010-09-22 


07:00:49 


1725 


378 


88.24 ±0.32 


NLTT17117 


2010-09-22 


08:42:38 


1800 


236 


88.61 ±0.40 


HIP74432 


2010-09-21 


23:16:51 


180 


211 


-39.03 ±0.40 


HIP74434 


2010-09-21 


23:23:45 


120 


120 


-38.56 ±0.36 


HIP81991 


2010-09-20 


23:31:06 


150 


270 


-5.29 ±0.25 


HIP94076 


2010-09-21 


23:32:49 


947 


439 


-44.08 ± 0.43 


HIP94075 


2010-09-21 


23:55:35 


2716 


548 


-40.08 ± 0.41 






UU:o4:U4 


on 

yu 


OiO 


— oo.oy ± u.zy 


HIP101932 


2010-09-21 


01:02:19 


1200 


451 


-52.43 ± 0.24 


HIP102532 


2010-09-22 


00:52:17 


30 


442 


-5.72 ± 0.32 


HIP102531 


2010-09-22 


00:59:06 


75 


443 


-6.56 ± 0.49 


HIP107759 


2010-09-22 


01:12:56 


3112 


286 


-63.95 ± 0.33 


HIP115126 


2010-09-21 


03:02:41 


190 


436 


8.49 ± 0.26 


NLTT56465 


2010-09-21 


03:16:26 


1478 


453 


11.55 ±0.19 


HIP114855 


2010-09-21 


01:52:46 


105 


624 


-25.62 ±0.28 


NLTT56278 


2010-09-21 


01:58:50 


2400 


276 


-22.88 ±0.36 
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Table 5. Measured Stellar Parameters 



Object 


Teff (K) 


logff [cgs] 


[Fe/H] 


Method^ 


hip9148 


5649 ± 27 


4.01 ± 0.03 


-0.146 ±0.020 


(1) 




5619 ±25 


3.93 ±0.06 


-0.165 ±0.022 


(2) 




5619 ±24 


3.94 ±0.05 


-0.164 ±0.022 


(3) 




5622 ±31 


3.88 ±0.06 


-0.172 ±0.024 


(4) 




5620 ±31 


3.94 ±0.05 


-0.164 ±0.022 


(5) 


hip9243 


5045 ± 17 


3.41 ±0.07 


-0.076 ±0.025 


(1) 




5048 ± 45 


3.46 ±0.06 


-0.068 ±0.025 


(4) 




5048 ± 45 


3.43 ±0.06 


-0.077 ±0.025 


(5) 


hip9247 


6135 ± 39 


4.12 ±0.06 


-0.193 ±0.025 


(1) 




6130 ±24 


4.20 ± 0.08 


-0.188 ±0.030 


(2) 




6130 ±24 


4.10 ±0.06 


-0.193 ±0.025 


(3) 




6131 ±41 


4.18 ±0.08 


-0.190 ±0.028 


(4) 




6140 ±41 


4.13 ±0.06 


-0.190 ±0.025 


(5) 


hipl0305 


6131 ±46 


4.13 ±0.06 


0.058 ±0.027 


(1) 




5996 ± 82 


3.80 ±0.05 


-0.022 ±0.041 


(2) 




5996 ± 82 


3.84 ±0.11 


-0.020 ± 0.039 


(3) 




6137 ±25 


3.86 ± 0.04 


0.035 ±0.057 


(4) 




6110 ±25 


4.09 ± 0.06 


0.046 ± 0.027 


(5) 


hipl0303 


5719 ± 14 


4.33 ±0.04 


0.112±0.015 


(1) 




5727 ±23 


4.31 ±0.18 


0.111 ±0.016 


(2) 




5727 ±23 


4.35 ± 0.05 


0.114±0.015 


(3) 




5739 ± 24 


4.31 ±0.18 


0.110 ±0.023 


(4) 




5742 ± 24 


4.40 ± 0.06 


0.125 ±0.017 


(5) 


hipl8824 


5880 ± 26 


4.07 ±0.05 


-0.066 ±0.027 


(1) 




5961 ± 108 


3.93 ± 0.06 


-0.042 ± 0.082 


(4) 
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Table 5 — Continued 



Object 


Teff (K) 


logs' [cgs] 


[Fe/H] 


Method^ 




5957 ±110 


4.25 ±0.07 


-0.025 ±0.030 


(5) 


hip20197 


4845 ± 26 


3.24 ± 0.06 


-0.035 ±0.028 


(1) 




4864 ± 49 


3.34 ±0.08 


-0.022 ±0.028 


(4) 




4865 ± 50 


3.30 ±0.07 


-0.031 ±0.027 


(5) 


hip23926 


5743 ± 24 


3.85 ±0.07 


-0.149 ±0.026 


(1) 




5703 ±61 


3.81 ±0.02 


-0.168 ±0.029 


(2) 




5702 ±61 


3.74 ±0.07 


-0.174 ±0.029 


(3) 




5706 ± 32 


3.81 ±0.02 


-0.168 ±0.028 


(4) 




5711 ±31 


3.77 ± 0.07 


-0.170 ±0.028 


(5) 


hip23923 


4963 ±110 


4.61 ± 0.03 


-0.171 ±0.079 


(4) 




4966 ±110 


4.69 ±0.14 


-0.130 ±0.085 


(5) 


hip32935 


4853 ± 16 


3.09 ±0.04 


-0.308 ±0.023 


(1) 




4916 ± 67 


3.23 ±0.08 


-0.281 ±0.024 


(4) 




4915 ± 67 


3.23 ±0.05 


-0.284 ±0.024 


(5) 


nlttiniT 


5835 ± 27 


4.24 ±0.07 


-0.405 ±0.030 


(1) 




5855 ± 51 


4.36 ± 0.07 


-0.382 ±0.033 


(2) 




5856 ±51 


4.29 ± 0.08 


-0.390 ±0.031 


(3) 




5772 ± 57 


4.32 ± 0.07 


-0.422 ± 0.047 


(4) 




5802 ± 58 


4.17 ±0.07 


-0.424 ±0.030 


(5) 


hip74432 


5686 ± 23 


4.13 ±0.05 


0.068 ±0.025 


(1) 




5690 ± 74 


4.28 ±0.05 


0.093 ± 0.040 


(2) 




5686 ± 75 


4.12 ±0.05 


0.074 ±0.025 


(3) 




5742 ± 49 


4.32 ± 0.04 


0.108 ± 0.034 


(4) 




5750 ± 49 


4.26 ±0.07 


0.099 ± 0.030 


(5) 


hip74434 


5573 ± 35 


4.30 ±0.10 


0.087 ±0.040 


(1) 
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Table 5 — Continued 



Object 


TeS (K) 


log 5 [cgs] 


[Fe/H] 


Method^ 




5573 ± 66 


4.11 ±0.13 


0.064 ±0.056 


(2) 




5577 ±64 


4.32 ±0.10 


0.084 ± 0.040 


(3) 




5737 ±86 


4.24 ±0.13 


0.106 ±0.090 


(4) 




5701 ± 84 


4.61 ±0.10 


0.137 ±0.053 


(5) 


hip81991 


5030 ± 16 


3.53 ±0.04 


-0.243 ±0.026 


(1) 




5061 ±64 


3.60 ±0.06 


-0.232 ±0.026 


(2) 




5061 ± 63 


3.60 ±0.04 


-0.235 ±0.026 


(3) 




5187 ±66 


3.66 ±0.04 


-0.223 ±0.057 


(4) 




5170 ± 64 


3.87 ±0.09 


-0.191 ±0.037 


(5) 


hip94076 


5953 ± 46 


4.33 ± 0.09 


0.151 ±0.031 


(1) 




5992 ± 20 


3.99 ± 0.08 


0.134 ±0.089 


(4) 




5951 ± 20 


4.33 ± 0.09 


0.152 ±0.031 


(5) 


hip94075 


5619 ± 23 


4.46 ± 0.05 


0.266 ±0.025 


(1) 




5489 ± 45 


4.40 ±0.10 


0.276 ± 0.054 


(2) 




5491 ± 45 


4.16 ±0.07 


0.214 ±0.040 


(3) 




5632 ± 18 


4.47 ±0.10 


0.272 ±0.025 


(4) 




5630 ± 18 


4.49 ± 0.05 


0.275 ±0.025 


(5) 


hipl01916 


5775 ± 14 


3.88 ±0.05 


0.068 ±0.016 


(1) 




5616 ±63 


3.68 ±0.02 


0.020 ±0.052 


(2) 




5611 ±64 


3.51 ±0.09 


-0.001 ±0.041 


(3) 




5660 ± 24 


3.71 ±0.01 


0.028 ±0.037 


(4) 




5666 ± 24 


3.62 ±0.06 


0.021 ±0.031 


(5) 


hipl01932 


5030 ± 33 


4.41 ± 0.06 


-0.007 ±0.023 


(1) 




5075 ± 19 


4.62 ± 0.02 


0.039 ± 0.043 


(2) 




5075 ± 20 


4.52 ± 0.04 


0.006 ± 0.027 


(3) 
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Table 5 — Continued 



Object 


TeS (K) 


log 5 [cgs] 


[Fe/H] 


Method^ 




5092 ± 33 


4.62 ± 0.02 


0.035 ± 0.039 


(4) 




5095 ± 33 


4.56 ± 0.04 


0.012 ±0.029 


(5) 


hipl02532 


4918 ± 33 


3.26 ±0.12 


0.142 ±0.056 


(1) 




4866 ± 47 


2.93 ±0.06 


0.061 ±0.074 


(4) 




4864 ± 46 


3.11 ±0.09 


0.133 ±0.055 


(5) 


hipl02531 


6403 ± 60 


3.91 ±0.07 


0.154 ±0.023 


(1) 




6277 ±67 


3.77 ±0.03 


0.088 ±0.039 


(2) 




6276 ± 67 


3.65 ± 0.08 


0.081 ±0.036 


(3) 




6297 ± 24 


3.77 ±0.03 


0.095 ±0.035 


(4) 




6301 ± 24 


3.69 ±0.08 


0.093 ± 0.032 


(5) 


hipl07759 


5530 ± 15 


4.03 ± 0.03 


-0.368 ±0.019 


(1) 




5437 ±60 


3.86 ±0.12 


-0.417 ±0.031 


(2) 




5437 ±60 


3.79 ± 0.08 


-0.427 ±0.031 


(3) 




5512 ±52 


3.84 ±0.12 


-0.403 ± 0.033 


(4) 




5495 ± 53 


3.97 ±0.05 


-0.390 ±0.021 


(5) 


hipll5126 


5557 ± 22 


4.13 ±0.05 


0.185 ±0.027 


(1) 




5383 ± 60 


3.86 ± 0.03 


0.142 ±0.047 


(2) 




5378 ± 60 


3.71 ±0.09 


0.110 ±0.046 


(3) 




5495 ± 30 


3.90 ±0.03 


0.135 ±0.038 


(4) 




5486 ± 29 


4.01 ±0.05 


0.153 ±0.028 


(5) 


nltt56465 


5113 ±56 


4.36 ±0.09 


0.181 ±0.036 


(1) 




5163 ±27 


4.56 ±0.03 


0.213 ±0.046 


(2) 




5163 ±27 


4.50 ±0.07 


0.192 ±0.039 


(3) 




5180 ±32 


4.57 ±0.03 


0.208 ± 0.044 


(4) 




5178 ± 32 


4.53 ± 0.06 


0.196 ± 0.040 


(5) 
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Table 5 — Continued 



Object 


Teff (K) 


logs [cgs] 


[Fe/H] 


Method 1 




hipl 14855 


4762 ± 38 
4766 ± 38 


2.76 ±0.04 
2.79 ±0.15 


0.109 ±0.060 
0.114± 0.060 


(4) 
(5) 


nltt56278 


4769 ± 90 
4772 ± 92 


4.58 ±0.02 
4.42 ± 0.16 


0.004 ±0.092 
-0.040 ± 0.082 


(4) 
(5) 



^(1) Excitation/ionization balance of iron lines, (2) Photometric 
Tcff and parallax, (3) Photometric TeS and ionization balance, (4) 
Ha Teff and parallax, (5) Ha T^g and ionization balance. 
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Table 6. Adopted Stellar Parameters 



Object Teff (K) logff [cgs]i [Fe/H] log^i^o [cgs]^ 



hip9148 


5625 


± 


35 


3 


.96 


± 


0, 


.07 


-0, 


.16 


± 


0, 


.04 


3, 


.90 


± 


0, 


.08 


hip9243 


5046 


± 


35 


3, 


.44 


± 


0, 


.06 


-0, 


,07 


± 


0, 


.04 


3, 


,49 


± 


0, 


.08 


hip9247 


6132 


± 


35 


4, 


.14 


± 


0, 


,07 


-0, 


,19 


± 


0, 


.04 


4, 


,15 


± 


0, 


.11 


hipl0305 


6109 


± 


69 


3, 


.93 


± 


0, 


,17 


0, 


,03 


± 


0, 


.05 


3, 


,77 


± 


0, 


.13 


hipl0303 


5730 


± 


35 


4, 


.35 


± 


0, 


,06 


0, 


,11 


± 


0, 


.04 


4, 


,30 


± 


0, 


.08 


hipl8824 


5893 


± 


68 


4, 


.06 


± 


0, 


,18 


-0, 


,05 


± 


0, 


.04 


3, 


,90 


± 


0, 


.07 


hip20197 


4854 


± 


36 


3, 


.29 


± 


0. 


.09 


-0. 


.03 


± 


0, 


.04 


3. 


.32 


± 


0, 


.09 


hip23926 


5716 


± 


36 


3 


.81 


± 


0. 


.05 


-0. 


.17 


± 


0, 


.04 


3. 


.88 


± 


0, 


.08 


hip23923 


4964 


± 


79 


4, 


.61 


± 


0. 


.05 


-0. 


.15 


± 


0, 


.09 


4. 


.62 


± 


0, 


.13 


hip32935 


4875 


± 


64 


3 


.16 


± 


0. 


.11 


-0. 


.29 


± 


0, 


.04 


3. 


.21 


± 


0, 


.14 


nlttl7117 


5828 


± 


53 


4, 


.28 


± 


0. 


.11 


-0. 


.40 


± 


0, 


.04 


4. 


.38 


± 


0, 


.13 


hip74432 


5711 


± 


55 


4, 


.23 


± 


0. 


.12 


0. 


.09 


± 


0, 


.04 


4. 


.30 


± 


0, 


.07 


hip74434 


5599 


± 


87 


4, 


.34 


± 


0. 


.23 


0. 


.09 


± 


0, 


.05 


4. 


.50 


± 


0, 


.07 


hip81991 


5077 


± 


92 


3 


.61 


± 


0. 


.11 


-0. 


.23 


± 


0, 


.04 


3. 


.57 


± 


0, 


.14 


hip94076 


5967 


± 


45 


4, 


.19 


± 


0. 


.26 


0. 


.15 


± 


0, 


.04 


3. 


.98 


± 


0, 


.13 


hip94075 


5587 


± 


86 


4, 


.42 


± 


0. 


.16 


0. 


.26 


± 


0, 


.04 


4. 


.29 


± 


0, 


.08 


hipl01916 


5685 


± 


83 


3 


.70 


± 


0, 


,11 


0, 


,03 


± 


0, 


.04 


3, 


,70 


± 


0, 


.07 


hipl01932 


5073 


± 


41 


4, 


.58 


± 


0, 


,09 


0, 


,01 


± 


0, 


.04 


4, 


,51 


± 


0, 


.08 


hipl02532 


4889 


± 


56 


3, 


.02 


± 


0, 


,19 


0, 


,12 


± 


0, 


.07 


2, 


,98 


± 


0, 


.13 


hipl02531 


6307 


± 


61 


3, 


.77 


± 


0, 


,08 


0, 


,11 


± 


0, 


.05 


3, 


,71 


± 


0, 


.08 


hipl07759 


5496 


± 


63 


3, 


.98 


± 


0, 


,11 


-0, 


,40 


± 


0, 


.04 


3, 


,84 


± 


0, 


.15 


hipll5126 


5488 


± 


83 


3, 


.93 


± 


0, 


,13 


0, 


,15 


± 


0, 


.04 


3, 


,90 


± 


0, 


.07 


nltt56465 


5165 


± 


37 


4, 


.54 


± 


0. 


.07 


0. 


.20 


± 


0, 


.04 


4. 


.50 


± 


0, 


.07 


hipll4855 


4764 


± 


35 


2, 


.76 


± 


0. 


.05 


0. 


.11 


± 


0, 


.05 


2. 


.71 


± 


0, 


.08 


nltt56278 


4770 


± 


66 


4, 


.57 


± 


0. 


.07 


-0. 


.02 


± 


0, 


.09 


4. 


.52 


± 


0, 


.08 
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Table 6 — Continued 



Object Teff (K) logg [cgs]i [Fe/H] log^iso [cgs]^ 



^Average surface gravity obtained as described in Section l3^ 

^Preferred value. Surface gravity obtained using T^s and [Fe/H] 
from this Table, observed visual magnitudes, Hipparcos parallaxes, 
and our isochrone technique described in Sectional 
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Table 7. Isochronal ages for wide binaries with turnoff and subgiant components 



ID 


"1" rl 1 G TXr/^vl^ 
tlilD WUl K. 






(Gyr) 


(Gyr) 


HIP 9148 


6.75ll;i 


7.776 ± 1.744 


NLTT 6583 






HIP 9243 


o nc + 0.67 
—0.59 


2.790 ± 0.423 


HIP 9247 


5.25^ 20 


5.594 ± 0.995 


HIP 10305 


9 -t.+0.60 
^- ' ^—0.33 


2.587 ± 0.326 


HIP 10303 


7 7c:+0.62 
' • ' ^—2.03 


4.832 ± 3.710 


HIP 18824 

NLTT 12412 


. -,^+0.47 


4.672 ± 0.483 


HIP 20197 
LSPM J0419+1419 


6.25t2-55 


5.559 it 1.573 


HIP 23926 

ilir 2:-!IJ2:-! 


- .-r 1 (i.')2 
"'■'"-2.34 


4.674 ± 0.332 

l.:-!Sl ± .--S.IJ.sO 


HIP 32935 


8.75l|°° 


7.891 ± 2.262 


NLTT 17117 




8.665 ± 3.655 


HIP 74432 


' ■ '3-1.78 


7.968 ± 2.045 


HIP 74434 


r, 9t-+4.92 
— 0.77 


8..384 ± 3.618 


HIP 81991 


A 7^+3. 56 
^- ' 3—0.06 


5.949 ± 1.636 


HIP 81988 






HIP 94076 


^•■^^—0.66 


3.490 ± 0.684 


HIP 94075 


°- '3-2.87 


4 AKQ -L Q Q9n 


HIP 101916 


o r,r; + 0.19 
O.ZO_(, 22 


2.944 ± 0.105 


HIP 101932 


6-25+3 J3 


3.655 ± 3.632 


HIP 102532 


1 9C+0.21 

^•^3_g 32 


1.114 ± 0.128 


HIP 102531 


9 7C+0.10 
^- '3-0.77 


1.774 ± 0.154 


HIP 107759 
LSPM J2149+1402N 




7.130 ± 2.977 


HIP 115126 




6.115 ± 0.869 


NLTT 56465 


3.25i5i^ 


3.699 ± 3.563 


HIP 114855 


1 9C+0.16 
^•^3_Q 34 


1.460 ± 0.260 


NLTT 56278 




5.021 ± 3.957 
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